Ultracold quantum dynamics: spin-polarized K + K 2 collisions 
with three identical bosons or fermions 
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We have developed a new potential energy surface for spin-polarized K( 2 S) + K2( 3 E+) collisions 
and carried out quantum dynamical calculations of vibrational quenching at low and ultralow col- 
lision energies for both bosons 39 K and 41 K and fermions 40 K. At collision energies above about 
0.1 mK the quenching rates are well described by a classical Langevin model, but at lower energies 
a fully quantal treatment is essential. We find that for the low initial vibrational state considered 
here (v = 1), the ultracold quenching rates are not substantially suppressed for fermionic atoms. 
For both bosons and fermions, vibrational quenching is much faster than elastic scattering in the 
ultralow-temperature regime. This contrasts with the situation found experimentally for molecules 
formed via Feshbach resonances in very high vibrational states. 
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I. INTRODUCTION 

Dilute gases of alkali-metal atoms are a rich source 
of novel physical phenomena. Bose-Einstein conden- 
sates (BECs) were first created in such gases in 1995 
[BSIjI and have been the subject of intense exploration 
ever since. Further possibilities were opened up by the 
achievement of Fermi degeneracy in 1999 Q. Among the 
alkali-metal atoms, lithium 0, Q and potassium 0, [|J 
have a special status because both bosonic and fermionic 
isotopes are available and both Bose-Einstein conden- 
sation and Fermi degeneracy have been achieved. The 
present paper will focus on potassium. 

In recent years, much interest has focussed on the in- 
teractions between atoms and the formation of molecules 
in ultracold gases. Donley et al. showed that it is pos- 
sible to form dimers of bosonic atoms such as 85 Rb by 
magnetic tuning from an atomic to a molecular state in 
the vicinity of a Feshbach resonance. However, molecules 
formed in this way in an atomic BEC proved to be short- 
lived (with lifetimes of milliseconds) because of atom- 
molecule and molecule-molecule collisions 0, H, H lioj : 
the molecules are formed in the highest vibrational state 
that exists in the two-body potential well, and any col- 
lision that changes the vibrational state releases enough 
energy to eject both collision partners from the trap. 

In summer 2003, Regal et al. succeeded in forming 
ultracold diatomic molecules in a Fermi-degenerate gas of 
40 K atoms by ramping the magnetic held through a Fes- 
hbach resonance. Such molecules are composite bosons. 
However, these too turned out to be short-lived (lifetime 
< 1 ms). Finally, at the end of 2003, a long-lived molecu- 
lar BEC was created using the same technique ^3 with 
a different Feshbach resonance. Such condensates have 
also been formed from dimers of 6 Li [l3l IhH ITU . 

The use of fermionic isotopes appears to be crucial 



for the production of long-lived molecular condensates. 
The inelastic collisions that cause trap loss for molecules 
formed from bosonic atoms are sometimes su ppr essed for 
molecules formed from fermionic atoms 0, Il7l fl8L fli^ . 
At magnetic fields where the atom-atom scattering length 
is large and positive (a > 1000 ao), molecular lifetimes 
longer than 100 ms can be achieved. Petrov et al. [2(j 
have explained the difference in collisional properties of 
dimers formed from bosonic and fermionic atoms in terms 
of the symmetries of the allowed wavefunctions. How- 
ever, their derivation applies only to dimers in Fesh- 
bach resonance states and not to deeply bound molecular 
states. 

In a previous study pl| . we investigated ultralow en- 
ergy collisions between spin-polarized Li atoms and Li2 
dimers, with the dimers in low-lying vibrational bound 
states. Our results showed no systematic differences 
in vibrational quenching between the bosonic 7 Li and 
fermionic 6 Li cases. This supports the conclusion that 
the suppression of inelastic collisions in the fermionic case 
requires both fermion symmetry |2(ij | and the long-range 
nature of the molecules in Feshbach resonance states [22|| . 

In the present paper, we study ultralow energy K 
+ K2 collisions involving three equivalent nuclei. We 
construct a new ab initio potential energy surface for 
the lowest spin-polarized electronic state of the potas- 
sium trimer (1 4 ^4 2 ) an d investigate both bosonic ( 39 K, 
41 K) and fermionic ( 40 K) cases. We perform quantum- 
mechanical scattering calculations at energies down to 
1 nK. Elastic, inelastic and rearrangement processes are 
considered. Our quantum dynamical results show that, 
for all three systems, vibrational relaxation is more ef- 
ficient than elastic scattering, as in our previous stud- 
ies with 23 Na [Hill and 6 ' 7 Li As in the case of 
lithium, there are no systematic differences in quench- 
ing rates for potassium dimers formed from bosonic and 
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fermionic atoms. 



II. POTENTIAL ENERGY SURFACE OF 

K 3 (1 4 A' 2 ) 

We have carried out ab initio calculations on K 3 
using a single-reference restricted o pen -shell variant 
|25| of the coupled-cluster method |2r| with single, 
double and non-iterative triple excitations [RCCSD(T)]. 
We used the small-core ECP10MWB effective core 
potential (ECP) of Leininger et al. [23 together with 
the medium-sized valence basis sets of Soldan et al. [28| . 
The quasirelativistic ECP treats the Is 2 electrons as 
core and the 3s 2 3p 6 4s electrons as valence. The valence 
basis set for K was used in uncontracted form. The 
resulting atomic polarizability (294.2 a n) is in excellent 
agreement with the experimental value [3(| (292.8 ± 6.1 

4). 

For interpolation purposes, the three-atom interaction 
potential was decomposed into a sum of additive and 
non-additive contributions, 

Vt I i mcI (r 1 2,r 2 3,r 13 ) = ^ Kumerfaj) + ^3(^2^23,^13)- 

i<j 

The full counterpoise correction of Boys and Bcrnardi 
|3l| was employed to compensate for basis set superposi- 
tion error in both dimer and trimer calculations. All the 
ab initio calculations were performed using the MOL- 
PRO package [13 • 

The dimer interaction energy Vdimer^) was first calcu- 
lated on an irregular grid of 42 points at interatomic dis- 
tances between 2.1 A and 14.0 A. The potential energy 
curve was generated using the modified ID reciprocal- 
power reproducing kernel Hilbert space (RP-RKHS) in- 
terpolation method [331 134| . The interpolation was done 
with respect to r 2 using RP-RKHS parameters m = 2 
and n — 3. Beyond 14 A the potential energy curve was 
thus extrapolated to the form 

v (*\ — —^1 — - o\ 

Vdimer(TJ - ^ 6 ^ 8 rW ■ (4) 

The long-range coefficients C§ and Cg were kept fixed 
to the published values of 3.897 x 10 3 E h a% and 4.2 x 
10 5 E h al respectively (HH^. The value of the "free" 
long-range coefficient Cio was then determined from the 
corresponding RP-RKHS coefficients [37j , and was found 
to be 2.0243 x 10 10 Eh aj°, which compares very well with 
2.0294 x 10 10 E h al° from Ref.H Values of D c = 252.6 
cm -1 and r c = 5.79 A calculated from the resulting curve 
are also in good agreement with experimental results on 
39 K 2 (D c = 252.74 ±0.12 cm" 1 , r c = 5.7725(20) A [H 

H). 

For quantum dynamics calculations, it is very impor- 
tant to have a potential energy function that can be in- 
terpolated smoothly (and without oscillations) between 
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FIG. 1: Cuts through the K3 quartet surface in valence coor- 
dinates. Upper panel: cut for a bond angle of 60°, showing 
the global minimum at —1269 cm -1 and 5.09 A. Lower panel: 
cut at collinear geometries; the collinear minimum is at —565 
cm -1 and 5.68 A. Contours are labelled in cm -1 . 



ab initio points. Oscillations often arise in low-energy 
regions if one or more points have much higher ener- 
gies than those surrounding them. Some experimenta- 
tion was needed to find a coordinate system in which 
interpolation could be carried out without problems. Ja- 
cobi and bond angle / bond length coordinates were 
rejected because they do not lend themselves to sym- 
metrization, and hyperspherical coordinates proved un- 
suitable because some combinations of grid points pro- 
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duce geometries with atoms very close together and cor- 
respondingly high energies. In the end, we chose to cal- 
culate the potential grid in pure bond length coordinates 
(^12) ?"23> f*3i)' This has the advantage that points that 
are related by symmetry have coordinates that are also 
simply related. 

The trimer interaction energy Vtrimer was calculated at 
325 points on a 3D grid covering the range of interatomic 
distances from 3.5 A to 10.5 A with step 0.5 A. To avoid 
duplication, only points with r 12 < r 13 < r 23 are re- 
quired. To meet geometrical constraints, all grid points 
must satisfy the triangular inequality \r\i — r\ 3 \ < r 23 < 
fi2 + fi3 . The distance r 2 3 was permitted to extend be- 
yond 10.5 A. The final grid consisted of 205 C 2v points, 
including 15 D3h points and 120 Coov points; the lat- 
ter include 15 Dooh points. Each calculation was carried 
out using the full symmetry allowed by MOLPRO. The 
non-additive energies V 3 were extracted from the trimer 
interaction energies using Eq. QJ. 



For low-energy scattering calculations it is important 
to have an interaction potential that behaves correctly 
at long range. The RP-RKHS interpolation procedure 
in 1 dimension allows this as decribed above. However, 
multidimensional RP-RKHS interpolation always gives a 
potential that extrapolates beyond the points as a simple 
product of inverse powers in the different coordinates. 
The leading long-range terms in the non-additive energy 
are the third-order dipole-dipole-dipole (DDD) jljj an d 
dipole-dipole-quadrupole (DDQ) terms given by 

T/ DDD _ o 7 (3) 1 + 3 COS 03 COS 01 COS 2 



v 3 DBQ = z{% (VP 123 + W 231 + W 312 ) (4) 

where 



TP 



ijk _ 



[9 cos < 



25 cos 30fc + 6 cos( 



4>j) x (3 + 5cos20 fc )] 



(5) 



and 4>i is the bond angle at atom i. It may be noted that the DDD term vanishes on a seam in the angular space 
and the DDQ term vanishes at all linear geometry configurations. Damped versions of these terms were therefore 
subtracted from the total non-additive energy V 3 before interpolation to give a quantity V3', 



Vl — V 3 — /damp [V 3 



DDD 



V, 



DDQi 



(6) 



The coefficients z[ 3 \ and Z^ 12 were taken to be 2.72 x 10 5 Eh a% and 5.11 x 10 6 Eh aj 1 respectively [42|. The damping 
function serves to prevent the non-additive energy exploding at short range, and was chosen to have a product form, 
/damp(ri2,f 23 ,ri 3 ) = f(r 12 ) /(r 23 ) f(r 13 ), where 



,(3) 



f(r) = exp - (k 3 /r - 1)' 



1 



< r < k 3 
r >k 3 



(7) 
(8) 



with the cut-off parameter £3 = 8.0 A. 

The leading term of the multipole asymptotic expansion of V 3 is the fourth-order dipole-dipole-dipole term (DDDD), 
which has a more complicated (unfactorizable) form [43j, 



V, 



DDDD 



45 



- 7( 3 ) 

- 64 Z iiii 



1 



1 



COS 2 <j>2 



1 



COS 



r 6 r 6 
' 12' 13 



12' 23 



~,6 r Q 
' 13' 23 



(9) 



(3) 

The coefficient Z\ 11X is not yet known, so this term cannot be subtracted out. However, the term is negative at all 

(10) 



geometries, so it can be eliminated by defining V 3 = jx V 3 , where 
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12' 23' 13 



(1 + cos 2 </>i) r^ 3 + (1 + cos 2 02 ) r? 3 + (1 + cos 2 3 ) rf 2 ' 

I 



The leading asymptotic term of the function V 3 now has 
the form —constant x t^z t y$ an( i ^ s suitable for an 
"isotropic" extrapolation of the type that results from a 



multidimensional RP-RKHS interpolation. The function 
V 3 was interpolated using the fully symmetrized 3D RP- 
RKHS interpolation method 44] . The interpolation was 
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v 3 = - y 3 " + U 



amp 



DDD 



DDQ 



(11) 



done with respect to the reduced coordinate {r/S) 3 and 
with parameters S — 10.0 A, m = 0, n — 2 in each in- 
teratomic distance. The original potential is then rebuilt 

as 

1 

9 

The final potential for quartet K3, Vt r i mer , has a global 
minimum at —1269 cm -1 at an equilateral (D3h) geom- 
etry 7*12 = T"i3 = ?'23 = 5.09 A. There is a shallow sec- 
ondary minimum at —565 cm" 1 at a linear Dooh geom- 
etry with T\2 = 7*13 = 5.68 A. Two cuts through the 
surface are shown as contour plots in Fig. ^ for values of 
the valence angle 60° and 180°. 

III. QUANTUM SCATTERING THEORY 
A. Method 

We have performed three-dimensional quantum dy- 
namical calculations for K + K 2 including reactive scat- 
tering for total angular momenta J = — 5. A time- 
independent formalism (which is the most appropriate 
choice for ultralow energy scattering) was used. The con- 
figuration space is divided into an inner and an outer 
region depending on the atom-diatom distance. In the 
inner region, typically for hyperradius smaller than p max 
= 60 ao, we use a formalism based on body- frame demo- 
cratic hyperspherical coordinates 0, which has pre- 
viously proved successful in describing atom-diatom in- 
sertion reactions such as N( 2 D) + Ho — » NH + H |47ll4c| 
and 0( 1 D) + H 2 -> OH + H gUHjJ. These coordinates 
were also used in our recent work on Na + Na 2 |23l Izif 
and Li + Li 2 [2~l| . 

At each hyperradius p, we determine a set of eigen- 
functions of a fixed-hyperradius reference hamiltonian 
Hq = T + V by expanding the wavefunction in a set of 
pseudohyperspherical harmonics. The reference hamilto- 
nian incorporates the kinetic energy T arising from defor- 
mation and rotation around the axis of least inertia and 
the potential energy V. A typical set of eigenvalues for 
K3 is shown in Fig. ??. At small hyperradius, the adia- 
batic states in each sector span a large fraction of configu- 
ration space and allow for atom exchange. The scattering 
wave function is expanded on this set of hyperspherical 
adiabatic states. This yields a set of close-coupling equa- 
tions, which are solved using the Johnson-Manolopoulos 
log-derivative propagator |51| . 

In the outer region, we use the standard Arthurs- 
Dalgarno formalism [53 based on Jacobi coordinates and 
assume that the atom-molecule interaction can be de- 
scribed by an isotropic potential U (R) . We thus assume 
that no inelastic scattering occurs in the outer region. 
The radial channel wavefunctions in this region are solu- 
tions of the equation (in atomic units) 

where p is the K-K 2 reduced mass and k and / are respec- 
tively the channel wavenumber and the orbital angular 



momentum quantum number. The potential U(R) be- 
haves asymptotically as — Ck-k 2 /^ 6 with a coefficient 
Ck-k 2 = 9050 EhdQ. Equation (|12[) is solved by a finite 
difference method in the range 60-10000 ao. The match- 
ing of the inner and outer wavefunctions is performed on 
a boundary which is an hypersphere of radius p = 60 ao. 
This yields the reactance K, scattering S and transition 
T matrices. 



B. Cross sections and Wigner laws 

At ultralow energies, for both bosonic and fermionic 
systems, only Ta elements with a relative atom-diatom 
orbital angular momentum / = contribute. We obtain 
elastic and quenching cross sections from the diagonal 
elements T„- of the transition T matrix, 



fc 2 



CT Q = -,(1-\1-T U \ 2 ). (13) 



The complex scattering length can be written as 



a = — lim ( 

2i fe^oV k 



(14) 



The elastic and quenching cross sections at threshold can 
therefore be written respectively in terms of the square 
modulus of the scattering length and the imaginary part 
of the scattering length, 



o E = 47r|a| 



4-7T 

-—Im(a) 
k 



(15) 



The elastic and quenching rate coefficients Ke and Kq 
are obtained by multiplying the cross sections by the 
atom-diatom relative velocity, k/p in atomic units, 



K E = k; 

A* 



K Q = -—Im{a). (16) 



At ultralow collision energy, where only / = is involved, 
the elastic rate coefficient vanishes as the collision energy 
decreases and the quenching rate coefficient becomes con- 
stant. 

The Wigner threshold laws j5^| for an orbital angu- 
lar momentum I give the following dependences on the 
collision energy for partial elastic and quenching cross 
sections, 



E. 



.11: 



E. 



1-1/2 

coll ■ 



(17) 



The energy-dependence for partial elastic and quenching 
rate coefficients is obtained by multiplying these cross 
sections by (2E co \\/ p) 1 ^ 2 , 



E 



21+1/2, 

coll ' 



K 



E. 



coll- 



(18) 



We also need state-to-state cross sections for the cal- 
culation of rotational distributions. For partial wave J, 
they are generated from the T matrix using the standard 
formula. 
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J+j J+j' 



a vj ^ v >j> - 2 ^(2 J + 1) l T <7'Z>jil 

W + ^j J=0 !=|J _ j|i , = | J _ i ,| 



(19) 



In this paper, we present for the first time differential cross sections (DCS) for an atom-diatom system at ultralow 
collision energy. The state-to-state magnetically averaged DCS is given by 



da 



vj—>v'j' 

dui 



l - y 

U + l) ^. 



4A&(2j + l) 



^(2J + \)d 3 m , . m (d C m)T^j, m ,. 



(20) 



where d J m , m is a Wigner reduced rotation matrix ele- 
ment. 

In Eq. H20fl . #cm is the scattering angle in the center- 
of-mass (CM) coordinate system. 6 cnl = 0° is defined 
as the direction of the CM velocity vector of initial K 
atoms and corresponds to forward scattering for the K 
products (and backward scattering for the K2 products). 
Backward scattering of the K products (forward scatter- 
ing for the K2 products) thus corresponds to 9 cm = 180°. 



C. Symmetries 

We are interested in elastic scattering and vibrational 
relaxation, which have different roles in the formation 
and stability of a molecular BEC. Elastic collisions are 
favorable for evaporative cooling towards condensation, 
whereas inelastic collisions provide a trap loss mecha- 
nism. For evaporative cooling, knowledge of the ratio of 
elastic to non-elastic rate coefficients is crucial. 

Here the system is composed of three indistinguishable 
atoms 39 K, 40 K or 41 K and the K3 potential energy sur- 
face is barrierless. We thus have to take into account 
two different collisional processes: elastic collision K + 
K2(w, j), and vibrational relaxation K + K2(v, j) — > K + 
K-2{v',j')- For this latter process, the rovibrational en- 
ergy E v iji of the product molecule is smaller than that 
of the initial molecule. 

We consider three atoms K in their stretched spin 
states, with F = F max = S + I and M F = F. For 
S = 1/2, this requires Ms = +1/2 and so the elec- 
tronic spin wavefunction is symmetric with respect to 
exchange of nuclei. 39 K and 41 K atoms have I = 3/2, so 
Mi = +3/2 in the spin-stretched state and the nuclear 
spin wavefunction is also symmetric. Since the total spin 
F is an integer, these atoms are composite bosons. For 
40 K, I — 4 so that F is half-integer and the atoms are 
composite fermions. For the spin-stretched state of 40 K, 
F = F max = 7 + 5 = 9/2 and M F = M max = +9/2. 
In this work, all magnetic interactions are neglected, so 
that the electronic and nuclear spin orientations are un- 
changed during the collision. 

The spatial part of the electronic wavefunction for the 
1 4 j4 2 state is antisymmetric with respect to exchange of 



nuclei. As a consequence, for 39 K or 41 K the spatial 
part of the nuclear wavefunction has to be symmetric, 
whereas for 40 K it has to be antisymmetric. In our quan- 
tum scattering code, this is achieved by selecting pseu- 
dohyperspherical harmonics in the basis sets for bosons 
and fermions to give the correct symmetry for the nu- 
clear wavefunction. The adiabatic states in each sector 
are obtained by a variational expansion on a basis of hy- 
perspherical harmonics with A± symmetry (of the com- 
plete nuclear permutation group S3) for bosons and A2 
symmetry for fermions. They are respectively fully sym- 
metric and fully antisymmetric with respect to particle 
permutations to account for the indistinguishability of 
atoms. 



D. Partial waves 

The behavior at the lowest collision energies is more 
subtle in the molecular case than in the atomic case. In 
atom-atom collisions, the fundamental difference between 
bosons and fermions is seen in the behaviour of cross sec- 
tions. For two identical spin-polarized bosons, the elas- 
tic cross section becomes constant as the collision energy 
goes to zero because only the s-wave (I = 0) contributes. 
In contrast, for two fermions, the elastic cross section 
decreases as the collision energy decreases because there 
is no s-wave and only the p-wave (1 = 1) contributes at 
ultralow energy. 

In atom-molecule collisions, by contrast, I = is al- 
lowed for both bosons and fermions. For three identical 
atoms, only one partial wave is involved at ultralow en- 
ergy: J = + for bosons and J = 1~ for fermions. These 
both include I = 0. 

The channels contributing to each partial wave J n up 
to J = 3 are shown for bosons in Table |T| The total 
mechanical angular momentum J = j + 1 is the vec- 
tor sum of the diatomic rotational angular momentum 
j and the orbital angular momentum 1. In this study, 
we consider initial 39 K2 or 41 K2 molecules initially in 
rovibrational state (v = 1, j = 0), so the parity II is 
{-iy +l = (-1)' = (-1) J . For 40 K 2 , only odd j is al- 
lowed so we consider the initial state (v = l,j = 1). 
The partial waves for fermions (with contributing chan- 
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TABLE I: Channels contributing to partial waves J , showing 
allowed values of the orbital angular momentum / for a rovi- 
brational state (v,j — 0) of 39 K2 or 41 K2 molecules formed 
from bosonic atoms. 



0+ 


1- 


2+ 


3" 


I = 


l = 1 


I = 2 


l = 3 



TABLE II: Channels contributing to partial waves J , show- 
ing allowed values of the orbital angular momentum I for a 
rovibrational state (v,j — 1) of a 40 K2 molecule formed from 
fermionic atoms. 



0+ 

I = 1 


1" 

/ = 0,2 


2+ 
1 = 1,3 


3" 
I = 2,4 




1+ 
I = 1 


T 
I = 2 


3+ 
I = 3 



nels shown in Table [H} can be separated into two cat- 
egories depending on the lowest value of I (or J7, which 
is the projection of the total angular momentum J onto 
the axis of least inertia in the corresponding hyperspher- 
ical basis set). The parity- favored partial waves include 
£7 = and the parity is the same as for the boson case. 
In contrast, for the parity-unfavored partial waves the 
£2 = component is forbidden (as is the lowest value of 
I) and the lowest value of £7 is 1. The parity- favored and 
unfavored partial waves have parity (— 1) J and (— 1) J+1 
respectively. 

E. Convergence and computer requirements 

K + K2 collisions are the most difficult and demanding 
ever studied in quantum dynamics calculations, because 
all three atoms are heavy and the potential well is deep. 
Careful attention has been devoted to the convergence of 
the calculations. There are three crucial parameters for 
convergence: the sizes of the basis sets, the size of each 
sector and the asymptotic matching distance. 

First, the number of fixed-p eigenstates included in 
solving the coupled equations in each sector must be 
large enough for convergence. It is essential to include 
many closed channels. For both bosonic and fermionic 
systems, this number, which is also the number of cou- 
pled equations, increases about from 250 for J = to 
1411 for J = 5. The time taken to solve the cou- 
pled equations varies respectively from 2 minutes to 4 
hours per collision energy on a Power4 P690 IBM com- 
puter. The fixed-p eigenstates are expanded in a pseudo- 
hyperspherical harmonic basis built from trigonometric 
functions, truncated at A max , the maximum value of the 
grand angular momentum. A max varies from 198 (867 
harmonics) at small hyperradius to 558 (6625 harmon- 



ics) at large hyperradius. The calculation to build the 
basis sets takes 180 hours and produces an output bi- 
nary file of 40 Gigabytes which contains all information 
for the close-coupling code. 

Secondly, the size of each sector must be small enough 
to give converged results. For both bosons and fermions, 
two different sector sizes were used: 0.025 ao from p = 8.0 
to 26.0 ao and 0.05 ao from 26 ao to the matching distance 
60 ao- This yields 1400 sectors, which is about 4 times 
larger than were used for the case of Na3 |23l |24| . 

Thirdly, the matching distance p max , where the adia- 
batic states arc projected onto the arrangement channels, 
must be chosen to describe correctly K2 molecules in the 
vibrational states required (in even j states for bosons 
and odd j states for fermions) . We have taken p max = 60 
ao, which is larger than the value of 50 ao for the Na3 
case. For bosons, the hyperspherical wavefunction was 
projected onto a set of 39 K2 or 41 K2 rovibrational func- 
tions with j max = (76,70,66,60,56,50,42,34,24,2) for 
v = 0, . . . , 9. For fermions 40 K2, rovibrational functions 
with j max = (77, 71, 67, 61, 55, 49, 43, 35, 25, 7) were used. 

At collision energies below 1 pK, only the partial wave 
that includes / = contributes to cross sections ( J = + 
for bosons and J = 1~ for fermions). However at higher 
energies, other partial waves contribute. In this work, 
for the bosonic cases we have converged the elastic and 
quenching cross sections for collision energies up to 10 
mK by including partial waves up to J = 5. Calculations 
for the fermionic system 40 K + 40 K2 are converged up 
to 100 pK by including partial waves up to J = 2. We 
have considered only three values of J for fermions be- 
cause the calculations are twice as expensive for fermions 
as for bosons (because both parity-favored and parity- 
unfavored partial waves contribute in the fermion case). 



IV. DYNAMICAL RESULTS AND DISCUSSION 

The following discussion refers to collisions involving 
39 K bosons unless otherwise stated. 



A. Hyperspherical adiabatic energies 

The hyperspherical adiabatic energies for the lowest 
250 hyperspherical states with £7 = are shown as a 
function of the hyperradius p in Fig. ?? (not shown here 
because the eps file is too big). If we consider 39 K2 
molecules in their (v = X,j = 0) rovibrational state, we 
have 15 hyperspherical states asymptotically open, and 
235 asymptotically closed. There are two minima: the 
deepest corresponds to the equilateral geometry, while 
the secondary minimum at around 20 ao corresponds 
to the linear geometry. At large p, the hyperspherical 
adiabatic energies tend to the energies of the diatomic 
molecule K2. The zero of energy is taken as the dissoci- 
ation limit to three atoms. 
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FIG. 2: Elastic cross sections for 39 K + 39 K 2 (t; = l,j = 0): 
partial as solid lines, total as dashed line. 
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FIG. 3: Quenching cross sections for 39 K + 39 K2(u = l,j = 
0): partial as solid lines, total as dashed line. Vertical lines 
correspond to the maxima of the effective potential (see text 
for detail). 



B. Cross sections and rate coefficients 

The fully converged elastic and quenching integral 
cross sections and their partial wave contributions for 
I = to 5 are shown in Figs. [21 and [21 for collision ener- 
gies between 1 nK and 10 mK. At E co \\ = 1CP 9 K, the 
elastic cross section is about 1.9 x 10 -12 cm 2 whereas the 
quenching cross section is about 1.4 x 10 -9 cm 2 . The 
slope of each partial cross section is given by the Wigner 
laws (Eq. i|17|) ) at ultralow energy except for the case of 
elastic scattering for I > 1, where the dispersion-modified 
threshold law a l E ~ £coii applies. The upper limit of the 
Wigner regime is 0.1 /iK. The minimum in the I = par- 
tial elastic cross section near 10~ 2 K arises simply from a 
near-zero phase-shift in the formula for the elastic cross 
section. 

The quenching cross section for a partial wave with / > 




10 9 10 



10 -7 1 6 10-5 10" 4 10- 3 

Collision energy (K) 



10" : 



FIG. 4: Total quenching rate coefficients for 39 K + 39 K 2 (v = 
l,j — 0). The result from the Langevin model is shown as a 
dotted line. 



shows a maximum at a collision energy which is given 
approximately by the maximum of the effective potential 
including the centrifugal and dispersion terms, 



V\R) = 



Ki+i) a 



The height of the barrier is 



K-K 2 



, [^ + 1)] 3/2 

maX 3^/2 ( 6Ck _k 2 ) 1/2 



at a distance 



R 



i 

max 



6^Ck-k 2 
1(1 + 1) 



1/4 



(21) 



(22) 



(23) 



The resulting barrier heights are included in Fig. |3J The 
first vertical line corresponds to the I — 1 partial wave 
and so on up to the / = 5 partial wave. It may be 
seen that each partial wave has a maximum at an energy 
slightly higher than the corresponding Vj' lax . At collision 
energies below the centrifugal barrier, the quenching par- 
tial cross sections for each I follow Wigner laws given by 
Eq. (|T7|) . Above the centrifugal barrier, the quenching 
probabilities come close to their maximum value of 1 and 
the cross sections vary as E~ x because of the k~ 2 factor 
in the expression for the cross section. 

The total quenching rate coefficient is larger than the 
elastic rate over a wide range of collision energies, up 
to 1 mK, and is three orders of magnitude larger than 
the elastic rate at 1 nK (Fig. [SJ . The large quenching 
rates are consistent with recent experiments that create 
molecules in atomic Bose-Einstein condensates by Fesh- 
bach resonance tuning 0, B IS HH) for bosonic atoms, 
such experiments were unable to produce a long-lived 
molecular condensate because of the efficiency of quench- 
ing collisions. 
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At high collision energy, when several partial waves 
are involved (I = — 5), the total quenching rate coef- 
ficient can be compared with that given by the classical 
Langevin capture model |54j|. 

KT^E) = *fe (24) 

This rate coefficient is shown as a function of collision 
energy in Fig. 0] It may be seen that there is semi- 
quantitative agreement between our quantum results and 
the capture model for collision energies above 0.1 mK. 
Below that energy, differences appear because there are 
fewer partial waves and the dynamics must then be de- 
scribed by a full quantum- mechanical treatment. 



C. Differential cross sections 
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FIG. 6: Differential cross section for vibrational relaxation at 
1(T 6 K (dotted line), 1CT 4 K (dashed line) and 10~ 2 K (solid 
line) . 
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FIG. 5: Differential cross section for elastic scattering at 
1(T 6 K (dotted line), 1CT 4 K (dashed line) and 1(T 2 K (solid 
line). The large forward scattering at 10 -2 K is not com- 
pletely shown; at zero degree the DCS value is 3 x 10~ 12 



Total differential cross sections (Eq. §'2Q\i ) f° r elastic 
and quenching processes with 39 K2 molecules initially in 
(v = l,j = 0) are shown in Figs. |S] and for collision 
energies I0" 6 K, I0~ 4 K and I0~ 2 K. At 10" 6 K, where 
only the I = and I partial waves contribute, the elastic 
and quenching DCS are quasi-isotropic. Below this en- 
ergy (not shown), only the / = partial wave is involved 
and so the DCS are fully isotropic. At higher energies, 
other partial waves contribute and the DCS depend on 
the center-of-mass scattering angle 9 cm . 

For the elastic process, Figure |S] shows that backward 
scattering is smaller at I0 -4 K than at 10 -6 K, whereas 
forward scattering is similar. At 10 -2 K, an enhancement 
of forward scattering is found. Undulations at this energy 
are simply due to a rainbow effect. 



In contrast, the behavior of the quenching DCS is 
less systematic (Figure EJ). The angular distribution is 
strongly peaked in the backward direction at 10 -2 K 
whereas the forward scattering is dominant at 10~ 4 K. In 
addition, at 10 -2 K the quenching DCS presents some os- 
cillations in the sideways scattering with two pronounced 
maxima at 30 and 70 degrees. 



D. Final rotational distributions 




FIG. 7: Rotational distributions as a function of the collision 
energy. The label j' is the final rotational quantum number 
of 39 K 2 (?/ = 0). 



The final rotational distribution for 39 K2 molecules ini- 
tially in (v = l,j = 0) is shown in Fig. as a func- 
tion of the collision energy. At each energy, we have 
divided each rotationally-resolved cross section (corre- 
sponding to a given final rotational quantum number j' 
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for 39 K2(v' = 0)) by the total quenching cross section. 
The sum of the distribution over j' gives unity at fixed 
energy. Because cross sections obey the Wigner laws (Eq. 
(|17|l . the rotational distribution becomes constant in the 
ultracold regime. The final rotational state j' = 8 is the 
most populated at 1 nK whereas j' = 2,4,6,8 are the 
most populated at 10 mK. Figure[7]also shows that only 
the lowest rotational states are significantly populated at 
all collision energies. 



E. Bosons versus fermions 

In Fig. [H] we compare the total elastic and quenching 
rate coefficients in the ultracold regime for the three sys- 
tems: 39 K + 39 K2 or 41 K + 41 K2, composed of bosonic 
atoms and 40 K + 40 K2, composed of fermionic atoms. 

In the ultracold regime, the total rate coefficients are 
very similar for the two bosonic systems. However the ro- 
tationally state-resolved cross sections (not shown here) 
are different and the similarity of the total cross sections 
and rate coefficients is evidently coincidental. The rota- 
tional distributions are sensitive to the mass and also to 
the interaction potential as in Na + Na2 |24fl . 

Figure |H1 shows that the quenching processes are very 
efficient even for the fermionic system. There is only 
a small difference between fermions and bosons. Both 
elastic and quenching rate coefficients are slightly smaller 
for the fermionic system than for the bosonic systems. 
However the ratio of elastic to quenching rate coefficients 
is nearly the same for bosons and fermions. 



10 mK, using a new quartet potential energy surface for 
the potassium trimer. We have found that the quench- 
ing rate coefficient is much larger than the elastic rate 
coefficient at ultralow collision energies. This applies 
for all three collisions studied here, 39 K + 39 K2, 41 K 
+ 41 K2, and 40 K + 40 K2 where the initial K2 molecules 
are in the (v — l,j — 0) rovibrational state for bosons 
or (v = l,j = 1) for fermions. Thus K2 molecules in 
these states are not good candidates to accumulate large 
densities of ultracold molecules or to achieve long-lived 
molecular BEC. For the bosonic systems ( 39 K and 41 K), 
the results are qualitatively similar to those we obtained 
previously for 23 Na [HEI and 7 Li [2l| atoms. 

The quantum results for quenching rates agree semi- 
quantitatively with the Langevin model at high collision 
energies, but below 0.1 mK the dynamics is described 
only by the quantum theory. 

The study of the fermionic system with 40 K atoms of 
particular interest, because it has recently been possi- 
ble to create stable molecular BECs of 6 Li2 and 40 K2 
fl2l ll3L IT1. Il5j . In these experiments, molecules were 
formed in the highest vibrational state supported by the 
potential well and it was found that inelastic collisions 
were suppressed by Pauli blocking [2(j ■ Our result shows 
that for 40 K2 molecules in the initial vibrational states 
v = 1, there is no suppression of the quenching process. 

Calculations with potassium atoms in higher vibra- 
tional states for both bosonic and fermionic systems are 
in progress. Future studies involving mixed collisions, 
i.e. with two different isotopes of potassium, are also 
planned. 
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FIG. 8: Comparison of elastic and quenching rate coefficients 
for bosons (dashed and dotted lines) and fermions (solid line) . 



V. CONCLUSION 



We have performed the first quantum dynamics calcu- 
lations for K + K2 collisions in the energy range 1 nK - 
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